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A Kubo inspired formalism is proposed to compute the longitudinal and transverse 
dynamical conductivities of an electron in a plane (or a gas of electrons at zero 
temperature) coupled to the potential vector of an external local magnetic field, 
with the additional coupling of the spin degree of freedom of the electron to the lo- 
cal magnetic field (Pauli Hamiltonian). As an example, the homogeneous magnetic 
field Hall conductivity is rederived. The case of the vortex at the origin is worked 
out in detail. This system happens to display a transverse Hall conductivity (P 
breaking effect) which is subleading in volume compared to the homogeneous field 
case, but diverging at small frequency like A perturbative analysis is pro- 

posed for the conductivity in the random magnetic impurity problem (Poissonian 
vortices in the plane). At first order in perturbation theory, the Hall conductivity 
displays oscillations close to the classical straight line conductivity of the mean 
magnetic field. 

PACS: 05.30.-d; 05.60.+w; ll.lO.-z 
1 Introduction : A classically free system 

Consider a Poissonian distribution of infinitely thin vortices carrying a flux 
(/), perpendicular to a plane (herafter called magnetic impurities). One would 
like to ask about the way such a system could be equivalent to a homoge- 
neous mean magnetic field, or, on the contrary, how it could localize electron 
wavefunctions, in such a way that their conductivity is altered. This jauestion 
might seem academic at first sight. However, it is commonly believedia that in 
Quantum Hall devices disorder does play a crucial role in the understanding of 
plateaus for the Hall conductivity as a function of 1/B or N{Ef), the number 
of electrons, at integer (or fractional) values in unit of e^/h. In the case of 
a homogeneous B field, linear response of the system to a small electric field 
gives no hint of such a remarkable behavior, since all states are delocalised and 
have the same transverse conductivity which varies linearly with 1/B or N[Ep) 
(classical straight line) . Disorder is needed to explain why some states (in fact 
most of them) are localised in broadened Landau levels, thus the plateaus in 
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the Hall conductivity, around the classical line. The random magnetic impu- 
rity model is a model where the disorder is contained in the definition of the 
magnetic field itself cl, and not added as an external random potential. After 
averaging over disorder, the model is translation invariant, as the usual Lan- 
dau system. Translation invariance is at the origin of the infinite degeneracy of 
Landau levels, which in turn is at the root of the understanding of the macro- 
scopic value of the Hall conductivity. We believe that the random magnetic 
impurity model encompasses at the same time generic features of a true Lan- 
dau system (and in a certain limit both systems are indeed equivalent) and of 
a disordered system. Finally, it is in essence a topological model, without any 
classical counterpart. But topological considerations □ have been put forward 
in the past to explain why, in the Quantum Hall Effect, although most of the 
states are in principle localised, still the quantized Hall conductivity remains 
macroscopic. 

It has been shown that two phases coexist in the random magnetic im- 
purity model: a phase where the average density of states of a test particle 
(an electron) coupled to the random magnetic impurity distribution exhibits 
Landau level oscillations, and a phase with no oscillation. We refer to the for- 
mer phase as an ordered phase, since the test particle "sees" a mean magnetic 
field irrespectively of disorder, and to the latter as a disordered phase, where 
the test particle "sees" impurities individually, and where disorder plays an 
important role. In the thermodynamic limit, the system is described by two 
parameters, a = where 0o is the flux quantum of the electron, and p, the 

mean impurity density. The natural interval of definition of a is [—1/2, 1/2], 
since the system is unchanged when the fiuxes are shifted by a multiple of the 
quantum of flux a ^ a + 1. In the absence of a physical orientation to the 
plane, as coming from an exterior magnetic fleld for example, reversing the sign 
of the flux a —>■ —a does not change the system either. So, one can restrict the 
interval of definition of a G [0, 1/2]. Note that the periodicity and the symme- 
try a ~a in the spectrum imply that the partition function is unchanged 
when a — > 1 — a, and therefore depends only on a(l — a). For a Poissonian 
distribution, with magnetic impurities randomly dropped on the plane without 
correlations -the simplest distribution considering the impurities as a Bose gas 
at zero temperature and introducing no additional parameter in the model-, 
one finds that the partition function of the mean magnetic field e{B) = hpa is 
reproduced as a series in power of (pa)", i.e. in terms of the average problem 
at n impurities. This result is a little surprising since it amounts to say that 
a purely free classical system -an electron coupled to Aharonov-Bohm fiuxes- 
leads, when quantized, to the quantum Landau problem. Strictly speaking, 
this limit is attained when a 0,p oo,pa fixed, which is singular. It is 



2 



possible to organize in a systematic way perturbative corrections to the lead- 
ing mean magnetic field expansion. It is sufficient to go at order p^a*, i.e. an 
electron interacting 4 times with 2 impurities, to get a critical value tto = 0.28 
below which one is certain that the average density of states has Landau level 
oscillations. 

What we are presently interested in is the transport properties of an elec- 
tron, or a gas of electrons at zero temperature, coupled to random magnetic 
impurities. We know for sure that in the mean magnetic field limit the con- 
ductivity is the Landau conductivity (classical straight line for the transverse 
conductivity) . We would like to see if the corrections due to disorder are going 
to produce oscillations -plateaus?- around the classical straight line. The pa- 
per is organized as follows: in the first section we will present a Kubo inspired 
formalisma where the dynamical linear response of a general 2d Pauli magnetic 
system to a small external electric field will be expressed in terms of the prop- 
agators of the unperturbed Hamiltonian. We will in particular emphasize that 
the time derivative of the dynamical conductivity is easier to compute than 
the conductivity itself, due to the particular structure of the Pauli Hamilto- 
nian. In the following section we will review, as a warm up exercise, the case 
of a homogeneous magnetic field (Landau problem), examine the difficulties 
inherent to the translation invariance of the system, and the zero width of the 
Landau levels, and present a way to circumvent them. Next, we will consider 
two examples of magnetic systems where the conductivity can be exactly com- 
puted : i) one electron or a gas of electrons coupled to one vortex, ii) one 
electron or a gas of electron coupled to a homogeneous magnetic field and one 
vortex. These systems are of particular interest since they both are related to 
the random magnetic impurity problem. In the next section, we will consider 
the random magnetic impurity problem itself, and use perturbative methods 
to get information on the Hall conductivity of this system. Finally, in the last 
section a discussion will follow where general conclusions will be drawn from 
the results obtained. 

2 Conductivity for Pauli magnetic systems 

By magnetic systems we mean in general the class of Hamiltonians for an elec- 
tron minimally coupled to a vector potential A{r) with the additional coupling 
of the electron spin up or down ct^ = ±1 to the local magnetic field B{r) (we 
set the electron mass = fi = 1) 

H^'-{p-eAir)y^^-^a. (1) 
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It rewrites 

<j, = +i i/„ = in_n+ (2) 
a, = -i i7d = in+n_ (3) 

where II-i- — (p^ — eA^) ± i(py — eAy) are the covariant momentum operators. 
In the homogeneous field case, the spin coupHng is a trivial constant shift, but, 
in general, it has important effects. In the one vortex or magnetic impurity 
cases, it is a sum_of 5{r — r^) functions, which is needed to define in a non 
ambiguous way □ H the short distance behaviour of the wavefunctions at the 
location of the impurities . 

One would like to evaluate the linear response of such systems to a small 
homogeneous external electric field in the x direction, i.e. to compute the 
longitudinal and transverse local conductivities Gxx {t, t) and Gyx (r, t) . Con- 
ductivities characterize the non equilibrium dynamics of the system under 
the influence of an electric field. Linear response theory relates them to the 
equilibrium correlation functions of the Hamiltonian H . In a time dependent 
formalism, one can start with an electric field E = 5{t)Eo, a single impulsion 
at initial time, to get a local current 

K^^\{Ar){r\ + \r){r\v} (4) 

{v is the velocity operator v = p — eA) proportional to the conductivity 

a,,ir, t) = ie{t)ei:T {p[j,{r, t),r, ]} (5) 

where 9{t) is the Heaviside function. Trjp- • •} is the thermal Boltzmann or 
Fermi-Dirac average. It can as well mean expectation value in a given quantum 
state if one wishes to compute the conductivity for an electron in this state, 
in which case p is the projector on this quantum state. ji{r,t) is the current 
density operator in the Heisenberg representation 

J(f,t) = e'^*J(f)e-'^* (6) 

We find convenient to compute conductivities in a propagator formalism. 
The propagator is defined as 

G0{f,r') = {r\e-''"\r') = ^ ^„(fX(f ') e"''^" (7) 

n 

where {If^Sn)} is a complete set of eigenstates for the Hamiltonian H. In 
the case of a thermal average for Boltzmann statistics, p is the usual density 
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operator p — c:xjp{—f3H)/Zf3, with 

Zp^ I drGf,{f,f) (8) 



the partition function. 

Consider the hnear combinations 

a^{'r,t) ^ axx{r,t) + i(7yx{f,t), {r,t) ^ axx{r,t) - iijyx{r,t) (9) 

The local thermal (Boltzmann) conductivity, (cr = cr^) 

ap{r, t)^Yl ^«(^' *)e"''''"/^/3 (10) 

n 

where a„{f,t) is the local conductivity for one electron in the state \(pn), 
rewrites in terms of propagators 

1 . 1_ 



'^^(^O = -ie{t)e'— j dr'{n±Gu{r,r)x'Gp-u{r',r) 



where IIj. is the hermitian conjugate of ll-i-. In (|Tl|) and in all formulae that 
follow, the differential operators 11 are always understood as acting on the r 
variable of the propagator that immediately follows them. 

The local conductivity (Jfj{T, t) leads upon space integration to the volume 
average global conductivity, for the global current density j — ev/V ^ 

a^{t) = ^e{t) ^TrMn±(t),x]} (12) 
In terms of propagators it rewrites as 

^(t) ^ ^ j dra^{r,t) (13) 

= ie{t)y^ J drdr'(u±Guir,r')x'Gf3-^tir',r)- {it ^ it + P) 

In the case of a system invariant by translation, the local conductivity cr/^^r, t) 
does not depend on r, and both terms in ( pj]) give the same contribution. 

To deduce from (|l^) the conductivity of a gas of electrons at zero temper- 
ature and Fermi energy Ep, one uses the integral representation of the step 
function 0{Ef - H) 

/oo ^iEpt' 
T^-T, —Zp^u'+e'(Tp^^t'+e'{t) (14) 
2z7r t' — irj' 



where e' and r]' are regulators which have to be set to zero at the end. The 
second-quantized conductivity rewrites as a sum of first quantized conductiv- 
ities because, as defined in (|^), it is a thermal average of a commutator of 
second quantized operators which are sums of first quantized operators. In the 
sequel we will refer to a^p as given in (|l^ as the "Fermi transform" of a p. 
One can rewrite m% as 



daE{t) ^ J_ 
dE ^ 2tt 

with 



dt'e^''''Zp^,t'+e'<yp^^t'+e'{t) (15) 



-^^(^) = 1/ '-^dE (16) 

In (p^, the integration starts at = since the spectrum is bounded from 
below at zero energy (Pauli Hamiltonian (|^)). 

In practice (Quantum Hall Effect devices), the number of electrons N{Ep) 
may be fixed, implicitly determining Ep by 



N{Ef) - TT- / dt'——^Zu'+e' (17) 



IT]' 

We also have for the density of states p{E) 



_ dN{E) 1 

with accordingly 



^(iS)s;:^ = ^/ <iiV«'Z„,+,. (18) 



In general, it will be more convenient to calculate the derivative &{t) of a{t) 
with respect to time, rather than a{t) itself. In the case of the thermal Boltz- 
mann conductivity, one gets 



a±(t) ^^l5{t)±9{t)^^ j drdr' (eB{T)n±Gu{r,r)x'Gp.,t{?' .^-iit ^ it+13)) 

(20) 



=2 

To derive (EG), the identity 



[i/, n±] = T |es(r)n± + [n±, ± y(f)] | (21) 



has been used, which is vahd in general for an Hamiltonian H = {p — eA)^/2 + 



V{r). But here, due to the particular structure of the Hamiltonian where 
V{r) — ~eB{r)az/2, we have 

[if„n+] = -eB{j^n+ (22) 

[Fd,n_] = eS(f)n_ (23) 



Therefore, ( pO| ) is valid if, when considering one computes , and, when 
considering H^, one computes . In both cases, we wish to emphasize the 
appearance of the local magnetic field i?(r) in (|2^). In the homogeneous field 
case, it factors out from the space integrals. In the magnetic impurity case, it 
is a sum of S{r — fi) functions, greatly simplifying the space integrals. 

We are looking at the Fourier transform of a(t) in frequency space, and 
more particularly at zero frequency uj — (i.e. the time independent or static 
response). The Fourier transform of a{t) is defined by 

/■oo 

CT(tj) = lim / dta(t)e^"^-'''>' (24) 

where e is a regulator which has to be set to zero at the end. The Fourier 
transform of ( pO| ) rewrites 

j drdf'{eB{r)Ii±Gu{ry)xG0.u{r\r)- {it ^ zt + /?))}(25) 
and an analogous formula for the conductivity of a gas of electrons cr^^ (uj) 

2 ( ^ pQC iEpt' pea p 

" (- - (-) = y [NiE,) ± ^ y_ ^t'y—^, I dt e(--)* y drdr' 

(^eB{r)U±Gu{r, r')x' _,t{r\ f) - (%t ^it + it' + e')) | (26) 

Up to now, the Hamiltonian (0) has been used. However, one can simplify 
the expressions above by noting that a potential vector A(r), in 2d, can be 
always rewritten up to a gauge as 

eA,(f) ^ ^,,d,m (27) 

with the covariant momentum 

H+ = -2i(aj + aj0) (28) 

H_ = -2i{d,,-dz<^) (29) 
7 



Let us redefine the wavefunctionsQ as 

Tpif) = U{f)-ip{r) (30) 

where the non unitary transformation U{r) respectively read in the spin up 
and down cases 

C/u(r) = e-*(^^) (31) 
Ud{r} = e+^^f^ (32) 

One has 

{/,-in+f/„ = (33) 

U-^H-Uu = U"_+Aid,<j) (34) 
U^'U+Ud = nl-Aid,(b (35) 
U^^n^Ud = Tl°_ (36) 

where the potential vector disappears from the covariant U,^^Il+Uu = Tl'^ 
and U^^n^Ud ~ n° operators, which narrow down to the free covariant 
momentum operators H'^ — —2idz and n"_ = ~2idz- The Hamiltonian H ~ 
U~^HU acting on Tp rewrites 

Hu = -2dzdz+^d(j>dz (37) 
Hd = -2dzd, - 4d,(t> dz (38) 

The propagator for the Hamiltonian H being defined in (Q) as Gp(f,r') = 
{r\e~^^\r'), one has 

GMr,f') EE {r\e-P"\r') ^ BL^G^ir,^ (39) 

It follows that the Hamiltonians iJ„, on the one hand, Hd, Hd, on the 
other hand, are equivalent, and can be indifferently used for computing ther- 
modynamical observables, i.e. traces of product of operators, as the partition 
function, the density of states, or the conductivity. The partition function 
rewrites as 

Z^^ I drG{r,r) (40) 



and the conductivity becomes 

a±(i) = ^e{t)^J drdr'{U-\m±U{7^G,t{r,r')x'GpMr',r) 

- {it it + P)] (41) 
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or 



V V Zip J \ 

- (it it + /3)) (42) 

It is remarkable that the covariant momentum operator U^^{f)IlU{r) appear- 
ing in cr or (7 in ( |4l] , ^ ) reduces to the free momentum provided that one 
considers cr+ in the spin up case, and a~ in the spin down case. It foUows 
that it is appropriate, using H, to compute (T+ in the spin up case, and in 
the spin down case. In the sequel, we will essentially concentrate on the spin 
down coupling. However, one should keep in mind that computations for the 
spin up case should follow the same lines. 

One finally extracts the real part of the transverse and longitudinal con- 
ductivity 

Re a^^iij) ± iRe a,,(^) = " ^"^^ ^'"^^ (43) 



(43) is a consequence of (^ and is valid both for the Boltzmann (|25|) or Fermi 
(26) cases. 

3 Hall conductivity for a homogeneous magnetic field 

As a warm up exercise, let us rederive the Hall conductivity a^p — 0) i'^ 
the homogeneous magnetic field case. By convention, and without any loss of 
generality, we assume eB > 0, and denote by ujc = +eB/2 half the cyclotron 
frequency. The Landau Hamiltonian is 

H.t - in^H^ ^^{p-LuJx f) ' - uj, (44) 
= Ih^H^ = ^{p-LuJxfy +u;c (45) 

where 11^ = -2i{d^- \ojcz) is the covariant Landau momentum. The constant 
spin couling shift ^lOc has clearly no influence on the conductivity. The shifted 
Landau propagator reads (-f(— ) corresponds to spin u (d)) 

Gg{r,r') = ^ce^'^'^" ^-iif(\z-z'\'' coth(puj,)+z'z-z'z) ^^q\ 
^ ' 27rsinh(/3cjc) 

Accordingly, the shifted partition function is 



^'^-^2fsinh(/5^c) ^^^^ 
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Its Fermi transform is the Landau density of states, namely in the spin up case 
p^{E) ^ - / dfe^^'z^,^, =V-Y1 - 2nc.,) (48) 

and in the spin down case 

I f°° 11°° 

p^{E) = — / dt'e^^' = E '^(^ - + (49) 

As announced above, let us concentrate on the spin down case, i.e. on 
a~ (with a homogeneous magnetic field, computing would lead in the spin 
up case to the same result). To compute the thermal Boltzmann conductivity 
(y), using 

n^G^(r,f') =iWc(^-^')(coth(/3wc) + l)G^(r,r') (50) 
one has to evaluate space integrals like 

I I dvdr\z - -z')G\,{r, f ')x'G^_,,(f', r) (51) 

Since, in the Landau case, Gff_{f,f')Gj^_^^{f',r) is a function of \z — z'\ only, 
( pT| ) might seem naively ill defined (infinite by power counting, vanishing by 
symmetry). Still, one can proceed if, looking at (pl|), one pays attention to 
first perform the r' integral, leading to the local conductivity a^{f,t), which 
upon integration on r, yields the global conductivity. In the homogeneous 
magnetic field case, the r ' integration in ( pl| ) gives, as it should, a result 
which is independent of r, the manifestation of the translation invariance of 
the system. The remaining r integral simply factorizes a volume factor, which 
gives an unambiguous meaning to the divergence in (|5^). One finds 

al^-it)=9it)ye^^'^^' (52) 



We would like to emphasize that (|52|) can be trivially recovered if one now 
looks at &j^~{t), instead of a^~{t). Indeed, in ( |20| ) B factors out from the 
space integrals. So, (|2^ rewrites as 

(t) =ySit) + 2iLo,a^- (t) (53) 

which can be directly integrated to (|5^), without any need of properly defining, 
in the thermodynamic limit, divergent space integrals. 
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In principle, one proceeds by taking the Fourier and Fermi transforms of 
(p^). The former reads 

^'i'^^^TF ./^o ^ ^{vV—^+n6i^ + 2u;A (54) 

V e — i(uj + 2ujc) V [ uj + 2ujc ) 



and the latter rewrites 



2 1 foo iEpt' 



4;W = ^W^e2-=*-l, / dt'^^-^Z,,>+,, (55) 
V 2m J i — irj' 

which yields the infinite series J^n^i^F — 2u!c{n + 1)). But this series is 
meaningless, because Ep is not defined due to the zero width of the Landau 
levels. 

However, coming back to (52) or (|5^), one sees that the thermal conduc- 
tivity does not depend on f3, implying that the conductivity is the same for 
each Landau state. It follows that, for a gas of electrons at zero temperature, 
the conductivity is obtained by multiplying (|2|) or ( p^ by N{Ef). From (p3|), 
one infers the transverse conductivity 

Re<(-)l,. = -iV(i^.)^^^ (56) 

In (|5^), the limit Wc ^ is properly defined only if one keeps w 7^ 0, in which 
case it vanishes, as it should. The Hall conductivity finally reads 

Re 4^ (c. = 0) I,, = ~N{Ep)^ 1 (57) 

This is the classical straight line, showing no plateaus in the Hall conductivity 
as a function of the number of electrons N{Ep), or of the inverse magnetic 
field 1/B. 

The longitudinal conductivity is 

2 

Re <y^^{co)Ux = N{Ep)y^ {d{u + 2lo,) + S{u; - 2w,)} (58) 



4 Hall conductivity for one vortex and for one vortex + a homo- 
geneous B field 

On the one hand, these two non trivial examples of magnetic systems are 
interesting because their conductivities can be entirely computed, thanks to 
the local magnetic field 6{-r) functions in the space integrals of o'/3(i). On the 
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other hand, they are respectively related, in the random magnetic impurity 
problem, to the conductivity of the mean magnetic field, and to perturbative 
corrections at first order in a to the conductivity of the mean magnetic field. 

As in the Landau case, we set eB > 0, thus ujc = +eB/2. We consider a 
vortex at the origin with flux 0. The coupling of the electron to the vortex is 
e(f)/2'!T — a. This system is periodic in a of period 1, therefore a can always 
be choosen in the interval [—1/2, 1/2]. If i? = 0, the interval of definition of a 
can be restricted to [0, 1/2], since reversing the sign of the flux does not change 
the system. If i? 7^ 0, can be either parallel (a > 0), or antiparallel {a < 0) 
to the magnetic field, which imposes an orientation to the plane. Clearly, the 
physics depends on the sign of a. We will focus on the physical situation where 
the flux and the magnetic field are parallel a G [0, 1/2]. Indeed, this is what 
happens in the random magnetic impurity problem, where the mean magnetic 
field (B) is built by the vortices (j) carried by the impurities. 

In the symmetric gauge, the Hamiltonian Hd reads 

— ' k yc v \ 

Hd = ^ [ p- uJck X f - a — 5— + {'Ka5{f) + uJc) (59) 



where k is the unit vector perpendicular to the plane. 

The coupling of the electron to the infinite magnetic field inside the vortex 
defines in a non ambiguous way the short distance behaviour of the electron 
wavefunctions at the origin. This need for a proper characterisation □ Q of 
boundary conditions at the location of the vortex is related to the fact that 



perturbation theory in a is ill defined for the Hamiltonian (59) in the zero 
orbital momentum m — sector. In fact, wavefunctions of zero orbital mo- 
mentum m = are affected by the short distance regularisation naSlf). In the 
case of the spin down Hamiltonian Hd, the contact term is repulsive, indicating 
that the wavefunctions have to vanish when r ^ as quickly as r". It follows 
that in the m = sector, only regular wavefunctions at the. origin should be 
retained, as in the standard Aharonov-Bohm prescription 0. It means that 
the electron cannot penetrate inside the vortex, a quite reasonable physical 
situation. 



4-1 Partition function and density of states 

The free and Landau partition functions diverge like the volume in the ther- 
modynamic limit (continuous spectra), and corrections due to the vortex are 
subleading in volume. To give an unambiguous meaning to these corrections, 
one needs to regularize infrared divergences. This can be achieved by adding 
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to the Hamiltonian an harmonic regulator, ■^LU'^r , which has to be set to zero 
tJo — > at the end. The Hamiltonian (59) with the additional harmonic term 
has the spectrum 

m ^ En,m = {'2n + \m — a\ + l)wf — (to — a)ujc + lOq (60) 

TO = -En.o = (2?^ + a + Ijti't + OLi^c + (61) 
where uj^ — uj'^ + lu'^. The partition function rewrites 

Z-°(i?,a)-Z-''(i3,0) = 



2smh{/3ujt) \ I - e-0('^t+'^'--') i _ e'^('^'-'-'-) 

(62) 

The thermodynamic limit, lOo — > 0, i.e. ujt — ujc ^ 0, yields 

«) - Z,iB, 0) = |. + ^3-^1 (63) 

The Fermi transform of ( |63| ) gives the density of states 

oo 

p{E,B,a) - p{E,B,0) ^ - ^(n + 1 - a) (5(£; - 2(n + l)w,) 

n=0 

oo 

+ ^ (n + 1) <5(£; - 2(n + 1 + a)iUc) (64) 

In the limit _B ^ 0, ( |63| ) becomes 

Z^(0,«)-Z;3(0,0)^ "(""^^ (65) 



and accordingly 



p(g,0,«)-p(i;,0,0)^ ^^ ^(g) (66) 



i.e. the usual Aharonov-Bohm depletion of statesEl at the bottom of the spec- 
trum with respect to the free density of states p{E, 0, 0) = p^{E). 

The physical meaning of ( ^4| ) is that on each Landau level En — 2(n + l)ajc, 
n + I ~ a states states disappear and n + 1 appear at energy 2(n + a + 1)lUc- 
When B = 0, the sole effect of the vortex in (|66| ) is that a(l — a)/2 states have 
left the bottom of the spectrum. 
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4-2 Hall conductivity for one vortex 

Let us first concentrate on the vortex alone. The propagator of ( p^ ) with B = 
reduces to the standard Aharonov-Bohm propagator 



where the /y(z)'s are the modified Bessel functions. 

One expects that the Hall conductivity is antisymmetric with respect to 
a —>■ I — a, because of the antisymmetry a — > —a (the transverse conductivity 
has to change its sign when the direction of the local magnetic field is reversed) 
and the periodicity. In particular it has to vanish when a — 1/2. One has to 
compute (pO|), where eB{r) = 2Tra5{f^ and n_ = —i{2dz ~ a/z). The 5{f^ 
appearing in cr'^{t) greatly simplifies a computation otherwise untractable. 
Indeed, only the m = and m — I orbital terms in the propagators eventually 
survive to yield 

^pit) = $m + $9it)-l-e^-'^'-^^^r it + ^py-"-t'-'^ {t-^Pr) (68) 
Its Fourier transform reads 

^^"(^) = Fr^-rll - ^e-^^^(r(l + a)n^ + «, 3; /3(c. + ^,)) 
^ 1/ (w + le) 1^ Zfj TT V 

-r(2 - a)«'(2 - a, 3; -/3(w + ie))) | (69) 

where the 4'(a,5;z)'s are the unregular confluent hypergeometric functions. 
To get the Hall conductivity Re (T(ti;)|j,2. in the w ^ limit (static response), 
special attention must be paid to the e — > 0'*' limit, because of the logarithm 
that appears in the low /3a; expansion of ^'(a, 6; I5{u} + ie)). One obtains (the 
electron mass, me, and h, have been reintroduced) 

„ , M ?ie^ sin(27rQ) / a{\ — ol) , „ ^^, ^ \ , s 

Re 'tM'^)!.. = ^2^2 ^2 [} + ^-^(.M^ ln(/3c.) + • • • ) (70) 

Thus, a single vortex in the plane is sufficient to display a Hall conductivity, an 
explicit P violating effect (seeB and references therein for early tentatives on the 
subject). Of course, it is subleading in volume a ~ XjV'^ when compared to the 
homogeneous _B- field case a ~ 1/1^- However, when, in the random magnetic 
impurity problem, comparaison will be made with the mean magnetic field, 
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the l/V"^ factor will be multiplied by the number of impurities N , yielding, in 
the thermodynamic limit, an adequate p/V behavior. Also worth mentionning 
is the leading 1/uP' divergent behaviour of (|70| ) in the — > limit, which will 
also be related to the dynamical conductivity of the mean magnetic field. 

For a gas of electrons coupled to the vortex at zero temperature, one has 
to compute 

At small frequency, one should look at the leading high t behaviour of (|68| ) in 
which /3 has been replaced by it' + e' . Since {t-t' + ie')^"" - t^~" {t + f - 
ie')" ~ —t', one finds a conductivity proportional to N{Ep) 

.,A^ - 0) . HE.fl [-^ + ^2.e-"sin(.a)^-l-^} (72) 



Thus, in the limit to ^ Ep, 

1 

,— sin(27ra)- 

and 



ReaEA^)\yx = N{Ef)tt^ siii{2TTa)— (73) 



ReaEA^)\xx = N{EF)^2sin^{TTa)^ (74) 



4-. 3 Hall conductivity for one vortex + a homogeneous magnetic field 
The propagator of the Hamiltonian ( |59| ) (spin down, a G [0, 1/2]) is 

27rsinh/3L^e ^ Uinh/3c^, 

(75) 

As in the one vortex case, and exactly for the same reason, only terms of 
angular momentum m = and m = 1 survive in {t) 

a-p{t) - -5{t)+2iu:,a-p{t) + -ie{t)^^—j—-^^G{t) (76) 
V ' V Z/3smh (pwc) TT 

where 

G(t) = d(t) {e**'^'^ (sinhiiwc)" (sinh(/3 - it)u;c)^"" - («* ^ + /?)} (77) 
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One can check that the hmit ^ of (t) in ( |76| ) gives (|68|), as it should. 
Its Fourier transform is 

. 1 e^r^ .1 cj^e-'5'^<=("+i) sin(7ra) ^ , ,1 

e - + 2wc) I smh (/3cjc) J 

where G(aj) is the Fourier transform of G{t). The thermal conductivity, splits 
in two terms, the usual Landau term for the magnetic field, and a term due 
to the vortex alone. Since G{t) is periodic, one expects the conductivity to be 
regular at low frequency. One has 

G(-) = YZ-J^ dt G{t)e-* (79) 

Since G{t)* — G{—t), G{lj) is purely imaginary. It follows that, in ([ts]), 
no contribution to the real part of the longitudinal conductivity has to be 
expected from the vortex. Furthermore, G{uj) is finite when w — > 0, since 
j^27r/we Q^^^ _ Q xhese two points are to be contrasted with the pure vortex 
case where the longitudinal and transverse conductivities have been shown to 
diverge at low frequency ( [73| , [74[ ). It means that, as in the pure Landau case, 
the limit cjc — > is not smooth when w — > 0. In other words, if one wants 
to recover ( |7^ , [7^ ) from the present analysis, one should, while keeping uj ^ 0, 
take the limit lUc —> 0, and then consider the small u expansion. Here, we are 
considering the other situation where lo —>■ 0, while Uc is kept finite. 

One is interested by the conductivity for a gas of electrons, i.e. the Fermi 
transform of ([78|), where /3 has to be replaced by a complex time (3 ^ it' + e' . 
It happens that the small a expansion of G{u! = 0) when (3 is imaginary has 
no term in a 

G{uj = 0) = i— smh{iujct') + 0{a^) + ... (80) 

Note, as a remark, that, when a — 0, G{lo) — i sinh(icL'ci')- 

It follows that the small a expansion of the Fermi transform of the static 
thermal conductivity (^) reads 



u; = = / <Zu' + -^{ h naHt' 



TTa 

2iTr t' - i-q' 2VuJc i"" ' V \ ~c 

(81) 

where the partition function Zn' given in ( |63|) should also be understood as 
expanded at first order in a. It is not surprising that in (^), the first order 
correction to the conductivity appears as a mean magnetic field contribution 
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of the vortex T^ajY to the magnetic field lOc 



(82) 

to be compared to the Landau conductivity (^3). 



5 Hall conductivity for the magnetic impurity problem 

5.1 The magnetic impurity problem 

Let us now onsider an electron coupled to the vector potential of a random 
distribution {fi, z = 1, 2, N} of N vortices at position r^, carrying a flux 0, 
hereafter called magnetic impurity. The thermodynamic limit is N —>■ oo,V —>■ 
oo, with a finite impurity density p = N/V. The distribution is Poissonian, i.e. 
point magnetic impurities dropped on the plane randomly without correlation. 
It means that the infinitesimal probability dPifi) of finding an impurity at 
position ri is 

dP{n) = ^ (83) 

One would like to see in which way this system could be compared with 
the naive mean field approximation where the local magnetic field eB{r) = 
27ra^^j^ 5{r — ri) is replaced by its mean value 

{B)V ~ Nip i.e. in the thermodynamic limit e{B) = 2tt pa (85) 

Two approaches have been used§ to study this problem: a path integral Brow- 
nian motion approach (or its discretised version, random walk in the plane, 
which allows for numerical simulations), which focuses on winding properties 
of Brownian paths in the plane, and a quantum mechanical Hamiltonian ap- 
proach, which relies mainly on perturbative expansions of the average partition 
hmction (small a, high temperature, etc). In the sequel, we will use the quan- 
tum mechanical Hamiltonian approach, that we summarize now. 



If one takes a 2 dimensional impurity density of the order of the density of current car- 
riers p = 4.10^^m~^, one obtains, for a = 1/2, a magnetic field precisely in the experimental 
range of the Quantum Hall Effect 

(B) = — ~ lOT (84) 
e 



17 



In the symmetric gauge, the Hamiltonian for N impurities is 



H^ = lip-Y. /P',f ) +-»T.^ir-n) (86) 



where a 6 [0, 1/2], without any loss of generahty. The need for a spin coupHng 
term has aheady been discussed above. Recall that, physically, it corresponds 
to hard-core impurities where the electron is unable to penetrate. 
Consider the nonunitary wavefunction redefinitions 



N 



iP{r^ = Y[\r-nr^ir) = Ud{r)^ir) (87) 



1=1 



become 



U-^n+Ud = Ul-2ian (88) 

u^'n^Ud = no (89) 

where the impurity vector potential ft = X^iLi ~ ^i)- The Hamiltonian H 
acting on ip^f) rewrites 

^ 1 

Hd = -29,5, -2a 9, (90) 

= ^n^ul- 2anii°_ (91) 

One can go a bit further by extracting the mean value of YiiLi \^ ~ ^il" 
which precisely yields, in the thermodynamic limit, the mean magnetic field 
Landau exponential factor e^^'^"'^'' with (wc) = e(i?}/2 — npa. Remember 
that a £ [0, 1/2], so that e{B), which has the sign of a, is indeed positive. So 
let us redefine 

N 



i=l 



and obtain 



{u'd)-'u+u'd = n<^> -2tain-{n)) (93) 
(f/^)-in_c/;i = n<^> (94) 
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Thus, 

H'a = ^n<-^>n<-^> - ia{n - {n))n<^> (95) 

where II^'^^ = —2i{dz — {uJc)z/2) is the covariant momentum operator in the 
symmetric gauge for the Landau Hamiltonian of the mean magnetic field, and 
(O) = Trpz is the mean value of f2. Therefore, (|95| ) is the Landau Hamiltonian 
for the mean field (B) up to a (ujc) shift, encoding the mean magnetic field 
contributions, plus a correction due to disorder, gathered in the interaction 
term V' = -ia{n - {n))Il<^> . 

The Hamiltonians Hd, Hd and H' d are equivalent, and can be indiffer- 
ently used for computing thermodynamical observables as the average partition 
function, the average density of states, or the average conductivity. However, 
quadratic interactions in the vector potential have disappeared from the Hamil- 



tonians 95), making easier the average over disorder. In the case of the 
Poisson probability distribution (^3|), the average can be done using identities 
like 



/ 



1 



z - Zi 



(96) 



dzidzi^^ ]_ _ = tt{^j^ + ^^) (97) 



What general information can be extracted from ( |8qj95| ) after averaging 
over the Poissonian disorder? 

i) because of the periodicity a — > a -f 1, and, in the absence of any ori- 
entation to the plane, because of the symmetry a — > —a, observables such as 
the partition function or the density of states, should be a function of a{l — a) 
only. The transverse conductivity, on the other hand, because of the electric 
field, should change its sign when a — > —a, implying that it vanishes when 
a = 1/2. 

ii) any observable can be expanded in power series in a and p as ^ p^a™" , m > 
n. The constraint m > n stems from the fact the number of vertex is neces- 
sarily larger than the number of impurities (an electron can interact at dif- 
ferent times with the same impurity). This series can always be rewritten as 

iii) the pa" and (pa)" terms are exactly known: the former correspond to 
the entirely solvable one impurity problem (standard Aharonov-Bohm prob- 
lem), whereas the latter is obtained, if, in the Hamiltonian iJ^, one sets 
fl = {fl), thus reducing the problem to the average Landau problem. In 
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the absence of disorder, any observable is necessarily a power series in leading 
order (pa)" of the mean magnetic field. 

iv) corrections coming from disorder are necessarily of the type p"'a"^ , m > 
n. To have information on these corrections, perturbation theory is needed. 
However, if perturbation theory is meaningful for the Hamiltonian (p6[), still 
ultraviolet divergences occur in the perturbative computationxiLthe spectrum, 
which are shown to cancel only if proper regulators are usedBB. On the con- 
trary, the Hamiltonians Hd or H'^, 95), yield a perturbation theory which is 



finite and easier to handle than whose of (p6|). If one concentrates on the per- 
turbative expansion for the Hamiltonian i?^, the unperturbed Hamiltonian is 
the Landau Hamiltonian for the mean magnetic field, that is to say (wc) = Trpa 
is fixed. It follows that p"a'", m > n terms correspond to a™"" corrections to 
the mean field. For instance, pa^ is actually of order a, meaning that a second 
order computation produces a first order a correction. Conversely, the a 
correction to the mean field will involve all the p"a"^^ terms. 

To develop a perturbative expansion for the partition function ( pO| ) or for 
the conductivity (41 4|), all what is needed is the perturbative expansion of 



the propagator G'p{r,r') — {r\e l^^d\r''^ of the Hamiltonian (|95|). At order 
n, one has 

5G'^'^\r, r') - (-1)" j^^ dPi j^' d/Js • • ■ 1^" ' d/3„ / dr i • ■ • / df„ 
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G<f > (f, r^)V' (r-i)G<i> ^ {r, , ri) • • • \/' {r^)G<^> {r^,?') (98) 



where V — —iaiVl — (i7))Hf and G^^^{f, r ') is the propagator of the mean 
magnetic field ( [46| ) for the spin down case with lUc ^ {'^c)- The only virtue of 
(r2) in V is to cancel the mean field (pa)" corrections coming from Q, alone. 
This is why one can ignore {VL) if one pays attention to keep only the non mean 
field contributions due to $7. In term of Feynman diagrams, it means that after 
averaging over the positions of impurities, one should only retain the diagrams 
without isolated impurity legs. It follows that the correction in a to the mean 
field expansion are obtained from the second order pa^ ~ {i-^c)Q' corrections, 
with diagrams involving only one impurity. 

Let us illustrate these considerations by computing the corrections to the 
average Landau partition function. One has to evaluate perturbatively the 
average second order correction to the average Landau propagator. ( p8| ) yields 
the correction 

SG'P (r, f') = {laf j^^ dp, J^^' J dridf^G<!:> (r, n) 
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N N 



Zi — z. 



Z2 - Zj 



X 



X 



(99) 



The average over impurity positions Vi and r}, using (|9q), produces either a 
term containing interaction with two different impurities i ^ j, a, p^a^ mean 
field contribution canceUcd by (fi) , or a term containing two interactions with 
the same impurity i — j, a pa^ term. It rewrites 



P 
X 



{6Gf\r,r')) 



r ^ 



(100) 



Zl- Z2 



= Trp{ia) / df3i / / dridr2 G p_o{f,fl)- 

Jo Jq J Z2 — Zl 

^-'^G<^%. (-^1 ' --'2)n<^> G<^> (f2 , r') 



Using the algebra of the average Landau operators 11! 



<L> ii<L> R<i> 



Bl 



and 



B<^>, with 



B 



<L> 



= -2i(di - -{ujc)z) 



B 



<L> 



= -2i{dz + -{uJc)z) 



(101) 
(102) 



and the completness relations for the propagators, the space integrations in 
(100) can be done to obtain 



{5Gf\r,r')) = -7rp(za)2^n<^>n<^>G<^>(r, f ') (103) 



(104) 
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Therefore, the correction to the average partition function at order a is 

FR) (/3(a;,))^ 

OZg = a 5 

' TT 2sinh2(/3K)) 

Accordingly, the correction to the average Landau density of states is 



(105) 



6p{E) ^a{LUc)- 



(106) 



n=0 



One can check that (105) reproduces the correction to the mean field term 



in the one vortex problem, by simply taking the small {lJc) limit in (105) 



(107) 



which is nothing but the term in ( |65| ) if one drops the pV = N factor. 
Considering the conductivity, we have in general to evaluate 



{it it + /?))^ (108) 



or its Fermi transform 



-2 roo ^iE^t' 



V 7„oo 2i7r t' - i-q' 



- / drdr'[li<'^>G'a{r,r')x'G'a'-u{r',r^ 



-it ^ t + t'))) 
At first order in perturbation theory, one has to compute 



(109) 



p x-^ 




+ 
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It is clearly a tedious task to evaluate perturbatively (108), already at first 
order in a. However, we will see below how to circumvent this computation. 



5.2 pa^ term: the one impurity case 

Before doing so, let us check that, at leading order in a, the conductivity 
for the one vortex case coincides with the low {lDc) expansion of the Landau 
conductivity for the mean magnetic field. This should be so, since we know 
that at leading order (pa)", the n random magnetic impurity system yields 
the Landau sytem ( [s^ ) at order (wc)" 

'.^X,^N{Ef)- — 



Rea<->{uo)\y,^N{EF)-^\l+[^ \ + ^ +••• (110) 



One has just to multiply the leading term in a of the one vortex conductivity 
( [73I ) by N to get the leading pa term 



ReaEA^)\vx ^ N{EF)^2^pa\=N{EF)\ye{B)\ (111) 



where the magnetic impurity density p — N /V, has factorized in (B). One 



recovers the leading term in the expansion (110) of the conductivity of the 
Landau problem in the low mean magnetic field limit (cjc)/w — > 0. If, in the 
single impurity case, the conductivity is divergent at small Pu), in the random 
magnetic impurity case, this divergence is nothing but a manifestation of the 
small {lOc)/lo expansion of the mean Landau conductivity. On the other hand, 
considering the longitudinal conductivity for the Landau problem (|58|), its 
small (wc) expansion starts at order (wc)f, a p^ behavior which cannot be 
obtained from the single vortex problem (|74|). 

5.3 p^-a"^^ terms: the one vortex + a homogeneous magnetic field 

As advocated above, perturbation theory with ( p5| ) at first order in a amounts 
to consider only the one impurity diagrams. It follows that, at this order, 
one can map the magnetic impurity perturbative problem on the exact one 
vortex + a homogeneous mean magnetic field problem. Consider again the 
Hamiltonian ( |5^ ) for the one vortex + (B) problem 

Hd = ^[p- {oJc)k X r - a'^r^ ] + {Tra5{f) + K)) (112) 
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and redefine the wavefunctions in order to extract only the short distance 
behaviour of the wavefunctions at the origin due to the vortex 



Ud{r) = r" (113) 

It follows that 

u;[^n+Ud = n<^>-2zai (ii4) 

u^^u^Ud = n<^> (115) 

The Hci Hamiltonian reads 

i/^^ in<^>nf^> -minf^> (ii6) 



Comparing the Hamiltonian (p3) with (116), one sees that they can be identi 



fied if, leaving asides the contribution of the average potential (Q), one restricts 
( p5| ) to one impurity located at the origin. This shows in particular that, apart 
for a factor N, the second order computation for the one vortex +{B) prob- 
lem is identical to the first order pa^ ~ (i-Oc)oi correction to the mean field in 
the magnetic impurity problem. One can check explicitly these considerations 



on the partition function. From (63), the term reads 



^ ' 2sinh2(/3(c.,)) ^ ^ 



Multiplying it by N reproduces the perturbative result (105) of the magnetic 
impurity problem at order a. 

5.4 Perturbative Hall conductivity 

We now go back to the magnetic impurity problem, and instead of computing 
perturbatively ( |108| ), we rather consider in ( ^l| ) the correction multiplied 
by N. One gets the conductivity for the magnetic impurity problem at first 
order in a 



2V{ujc) { ' ' ' ' dEp 

where Z^^-^ and N'^^^ (Ep) are respectively the partition function and the 
number of electrons for the average Landau problem. N'^-'^^ (Ef) is nothing 
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but the mean field contribution to the actual number of electrons in the system 
N{Ef) = N<^>{Ef) + aN^^^Ep) + a^N^^^Ep) + ■■■ where the iV(")'s are 
functions of (lUc) and Ep only. Then, one can rewrite 



<Jp^{oJ = 0) = ^^^N{Ep + a{^,)) (120) 

which is still valid at first order in a. The Hall conductivity is thus, at this 
order, 

Re apA^ = 0)1,, = -N{Ep + a(^,))l^ (121) 

a quite simple result (eventhough higher order perturbation theory might 
change this situation). 

Also, from (|7^), one finds no correction for the longitudinal conductivity, 
which means that at this order, ( p8| ) is not affected by disorder. 

6 Heuristic considerations and open questions 

Let us consider the experimental situation where the magnetic field (B) is kept 
fixed, and where the electron density varies. The HalljXonductivity ( |121| ) is, 
basically, N{Ep + a{ujc)) as a function of N{Ep). In El, the average density 
of states p{E) was estimated on analytical (approximation) and numerical 
grounds. In Figure |l|, p{E) is given for a = 0.01. This curve has essentially 
an a behavior, with broadened Landau levels having a width of order a{LUc), 
a height of order 1/a, such that the total number of quantum states per unit 
surface in a broadened Landau level, (lUc) /tt, is conserved whatever the disorder 

is. 

In Figure H, the Hall conductivity (121) is displayed (full line) as a function 



of the number of electrons N{Ep), while keeping the mean magnetic field fixed. 
The Hall conductivity has oscillations close to the classical straight line of the 
average Landau conductivity. We believe that the oscillations are smoothed 



because of the approximation made in (12C) and in p{E), and consequently 
in N{Ep). In fact p{E) is too broad to get a stepwise Hall conductivity. On 
the other hand, we also considered that p{E) at small a is the same object as 
p{E) at first order in a. But the expansion of the density of states around the 
mean field density of states is surely not analytical. To conclude, the curve 
given in Figure ^ should only be viewed as a crude first step towards a more 
complete computation. It is not impossible that higher order corrections, or a 
more general argument, will yieid a behavior closer to the experimental curve 
(dashed lines in Figure ^, fromB). 
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What kind of physical picture could support the model presented above? 
When the magnetic field is very strong, a ~ 1/2, the system is totally disor- 
dered, and cannot possibly transport any current. This is consistent with the 
symmetry argument which indicates that at a = 1/2 the transverse conduc- 
tivity has to vanish. On the other hand, when (B) decreases, Landau level 
oscillations appear, implying that some states conduct, thus the appearance of 
plateaus. When (B) becomes smaller and smaller, i.e. a —>■ 0, a, pure Landau 
system is reached with a classical Hall conductivity and no plateau at all. 

Clearly, no edge or finite size effects have been considered whatsoever, 
and no interaction between electrons either. Also, the thermal average has 
forbidden any precise information on the possible localisation of quantum states 
in the random magnetic impurity distribution. 

It would be certainly rewarding to know how to extract from (108) a 
more global (non perturbative) information, which would, like the Thouless 
topological argument □, tell us for sure that, when considering the Fermi gas at 
zero temperature and Fermi energy Ep, the correction terms always oscillate 
around the classical straight line. 
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Figure 1: Average density of states of the random magnetic impurity model for a = 0.01 
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Figure 2: Hall conductivity in unit oi /h of the random magnetic impurity model at first 
order in a for a = 0.01 as a function of the filling factor u = ^(f p"> ^ ■ straight line = 
classical result, steps = experimental Integer Quantum Hall Effect (schematic), full line = 

perturbative result 



28 



